HIP-2002-26/TH 



(N 
O 
O 
(N 

O 

Q 



> : 
o . 
m ■ 

CN ■ 

: 

o . 

CN ■ 

o ; 
^ : 

Oh! 

i 

CU. 
Q . 
^ ■ 



X 



Dependence of hadron spectra on 
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Using equilibrium hydrodynamics with initial conditions for the energy and net baryon number 
densities from the perturbative QCD + saturation model, a good simultaneous description of the 
measured pion, kaon and (anti)proton spectra in central Au+Au collisions at yfs = 130 AGeV is 
found with a single decoupling temperature Td cc = 150 . . . 160 MeV. The interplay between the 
resonance content of the EoS and the development of the transverse flow leads to inverse slopes and 
(Pt) of hadrons which increase with decreasing Tdcc- The origin of this result is discussed. 

PACS numbers: 25.75.-q, 12.38.Mh, 47.75+f, 24.85.+p 



Introduction. The amount of initially produced matter 
in ultrarelativistic heavy ion collisions at collider energies 
has been suggested to be controlled by gluon saturation 
|Q, ||, ||. For heavy nuclei, A ~ 200, the saturation scale 
which determines the dominant transverse momentum 
scale is expected to be 1 ... 2 GeV, which can lead to 
quite large values for the initial transverse energy [g, ||, |) . 
For instance, the initial dEx/dy at y ~ obtained from 
the perturbative QCD (pQCD) + saturation model Q 
exceeds the value observed in central Au+Au collisions 
at Brookhaven RHIC by a factor 2.5 ... 3 [g g. 

A mechanism to transfer transverse energy into the 
longitudinal motion is provided by the asymmetry of 
the initial collective motion of produced matter. The 
observed rapidity distributions of final particles suggest 
that the matter is produced in a state of rapid longi- 
tudinal expansion Q. In the collective expansion en- 
ergy is transferred through the work by pressure, pdV. 
Hence the strong initial expansion in longitudinal direc- 
tion leads to a large transfer of energy from the transverse 
into the longitudinal motion. The space-time evolution 
of a dense matter with collective effects and with a QCD 
phase transition is describable in terms of relativistic hy- 
drodynamics § § § @. 

Hadron spectra are often argued to be insensitive to 
the early stages of the collision. This is not quite true for 
events with large multiplicities since secondary collisions 
will lead to transverse collective motion and correlations 
in the properties of particle spectra which are difficult to 
understand otherwise. A clear correlation resulting from 
radial flow is the mass dependence of slopes of transverse 
spectra. For non-central collisions hydrodynamics pre- 
dicts the experimentally observed elliptic flow pl[ . 

As energy is transferred into longitudinal direction 
during the expansion, it is usually expected that the 
transverse spectra of final pions become steeper at lower 
decoupling temperatures. However, since the number 
of hadrons and resonance states increases rapidly with 
mass, a considerable fraction of the energy of thermal 
matter can be in the form of heavy resonances. When 



the total transverse energy of all hadrons decreases with 
the decreasing decoupling temperature, the energy re- 
leased from the reduction of the number of resonances 
and heavy particles can lead to an increase of the slope 
of the spectra of the remaining final stable hadrons. The 
energy release from the latent heat in a phase transition 
has a similar but smaller effect. In this work we study the 
details of these effects and the interplay between them by 
comparing the expansion of a massive pion gas with that 
of a hadron resonance gas. 

Slopes of hadron spectra from heavy ion experiments 
at Brookhaven AGS and CERN SPS, can be reproduced 
with a kinetic freeze-out temperature ~ 120 MeV [ p"2| . 
To reproduce the strange particle abundances in thermal 
models, chemical freeze-out temperatures higher than 
120 MeV are needed Q. Also at RHIC, a higher de- 
coupling temperature is supported by the heavy particle 
yields Q. In this note we argue that in central Au+Au 
collisions at yfs = 130 ^4GeV also the kinetic freeze-out 
takes place effectively at higher temperatures: solving 
the transverse expansion within boost invariant hydro- 
dynamics with initial conditions for the energy and net 
baryon number densities from the pQCD + saturation 
model [^, we describe simultaneously the measured 
7T^~ , and p(p) spectra with a single decoupling tem- 
perature Td e c = 150 . . . 160 MeV. Similar observation has 
been made in Ref. |14|j using a fireball parametrization. 

We describe briefly our framework to calculate the 
hadron spectra and then discuss how the hadron and 
resonance content in the Equation of State (EoS), and 
the phase transition affect the dependence of the slope of 
the spectra on the decoupling temperature. Finally, the 
pion, kaon and (anti)proton spectra are compared with 
the experimental results from PHENIX @ [ij, 0. 

The theoretical framework. The calculational frame 
which we use has been discussed in detail in The 
initial particle production is calculated from the pQCD 
+ saturation model based on lowest order pQCD interac- 
tions and a cut-off scale p S3 ,t determined from a saturation 
condition for the final state minijets m. The iVfactor 2.3 
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has been fixed on the basis of the next-to-leading order 
calculation of minijet transverse energy Q. 

The transverse energy weighted minijet cross section 
in a central rapidity unit, o{Et) (p S at, A, \y\ < 0.5), is 
the key quantity for determining the initial energy den- 
sity at the time of formation of the matter. As a new 
ingredient, we include the corresponding quantity for the 
initial net baryon number, <j(Nb) = (a(N q ) — a(Nq))/3, 
where the computation of a(N q ^} is based on the flavour 
decomposition of the minijet cross sections in |19[ |. 

Below, we shall consider Au+Au collisions at y/s = 
130 AGeV with a 5 % centrality cut. As explained in 
this corresponds to central collisions of effective nuclei 
with A e g = 181. The number of participants which we 
get for such collisions is 346, consistent with the PHENIX 
result 348 ± 10 Jl5fl . For the computation of the initial 
densities, we obtain <j(Et) — 67.0 mbGeV, cf(Nb) = 
0.64 mb and p sat = 1.06 GeV. Thus, the initial E T = 
1750 GeV and N B = 16.7 in the interval \y\ < 0.5 Q. 

The matter is assumed to be (approximately) thermal- 
ized at the time To = l/p S at = 0.19 fm/c. For cen- 
tral collisions, this approach has predicted successfully 
H the particle multiplicities at RHIC Q and also, 
when amended with a hydrodynamical description of the 
transverse expansion of produced matter, the transverse 
energies of final particles, dE T /dy, within 10. . . 20 % @. 

The high temperature phase of our EoS is QGP with a 
bag constant B and the low temperature phase a hadron 
resonance gas (HRG) with hadrons and hadron reso- 
nances up to M = 2 GeV and a repulsive mean field 
among the hadrons The values of the bag constant 
B and the mean field strength K are so chosen that the 
transition temperature is T c = 165 MeV. 

The thermal spectra of final stable hadrons are ob- 
tained by folding the thermal motion with the fluid mo- 
tion using the Cooper and Frye prescription g, pl| . 
The dependence of the spectra on the transverse flow 
velocity v r = tanhy r and on the decoupling tempera- 
ture Tdcc enters through the factors (pr/Tdec) coshy r and 
(fBT/Idec) sinhy r as the arguments of modified Bessel 
functions E.g., for pions with pt^I GeV, con- 

tributions from regions where i> r >0.5 behave approxi- 
mately as ~ exp(-p T /e yr T dcc ) = exp(-p T /T cff ) with 
Toff = e Vr Td cc . This shows explicitly how the inverse 
slope, T e ff, depends at large pr on the decoupling tem- 
perature and on the transverse flow rapidity y r . Due to 
the exponential dependence on y r , a small increase in the 
flow can compensate the decrease in Td ec - The net effect 
on T e ff , as the decoupling temperature and the transverse 
flow change, depends on the EoS and, in particular, on 
the assumed particle content of the hadron phase. 

To obtain the spectra of final stable hadrons, we first 
calculate the spectra of all hadrons and hadron reso- 
nances which are included in the EoS. Then the decay 
contributions from all resonance states are added to the 
spectra of stable hadrons. All strong and electromagnetic 



two and three body decay modes with known branching 
ratios are included ]22] |. 

dN(y=0)/dydp x 2 [GeV" 2 ] y r - y u 




Pt [GeV] lo g (T,/T) 



FIG. 1: Left: Transverse spectra of thermal n for two de- 
coupling temperatures T^ cc in the four cases studied; decay 
pions are not included for (b) and (d). The spectra for (c) and 
(d) have been scaled. Right: The change of the transverse 
flow rapidity y r as a function of temperature along a flow 
line starting at r = 5 fin at tq. The flow rapidities y Tt i at 
7i = 160 MeV are indicated. On the dashed-dotted line the 
effective temperature T c g = e Vr Tdc C would remain constant. 

The effect of resonances and phase transition on Td ec 
dependence of spectra. Before comparing our results with 
data, we try to clarify the Tdcc dependence of final hadron 
spectra by comparing the spectra of thermal pions at 
Tdcc = 160 and 100 MeV in four different cases: (a) mas- 
sive pion gas (PG) without phase transition (PhT), (b) 
hadron resonance gas (HRG) without PhT and (c) PG 
with PhT to QGP, (d) HRG with PhT to QGP. To make 
the differences in the effects on flow more transparent, we 
first consider only thermal pions also in cases (b) and (d) . 
The spectra of positive pions are shown in the left panel 
of Fig. 1 for each case. Except in the case (a) the slopes 
in the interval 1-5<pt^4 GeV change considerably. Pi- 
ons in this part of the spectrum come mainly from matter 
with large transverse velocity where T e ff = e y ' Tdec holds. 

In the expansion the temperature of the matter drops 
and the transverse velocity increases. These changes have 
opposite effects on the slope. They nearly cancel if the 
growth of the flow in terms of an increase in y r satisfies 
&Vr = Ur — yr,i = l°g(Ti/T) when the temperature drops 
from T\ to T . In the right panel of Fig. 1 we plot Ay r 
as a function of log(Ti/T) between T x = 160 MeV and 
T = 100 MeV for a fluid element moving along a flow 
line which initially starts at r — 5 fm. This fluid element 
belongs to the region which dominates the tail of the 
spectra. We see that for PG without PhT, the case (a) 
with the small change of slope (cf. Fig. 1, left), the curve 
is close to the dashed-dotted line for e Vr Ta cc = const. 
For the other cases the increase in flow with the decrease 
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111 


0.54 
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0.55 
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0.44 
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570 
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204 


0.96 
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0.57 
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0.56 
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0.54 
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0.51 
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588 
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319 


0.61 
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0.61 
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272 


0.56 
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TABLE I: Multiplicity dN* /dy and average transverse momentum (pf. }/GeV of positive pions for two different decoupling 
temperatures T e ff/MeV and the same four EoS as in Fig. 1. Resonance decays (RD) have been added in the case (e). The 
multiplicity dN/dy and the transverse energies {dEr / dy) / 'GeV and (dEr / drj) / GeV for all particles are shown for the case (e). 



in temperature is stronger and the changes in Ay are 
comparable. This indicates a similar relative increase in 
T e g, the inverse slope of the spectra, in the cases (b)-(d). 

The development of transverse flow is fastest for PG 
with large pressure gradients and slower for the other 
cases with softer EoS. In the cases (b) and (d) with HRG, 
a considerable amount of energy is stored in the heavy 
hadron and resonance states. As the matter expands, 
the release of energy from these states slows down the 
decrease in temperature leading to a stronger growth of 
the flow, as a function of decreasing temperature, than 
in the PG of the case (a). In the case (c), the latent 
heat released at the strong first order PhT leads to a 
similar behavior of the y r {T) as in the case (b) and (d). 
The reason why the change of slopes between (b) and (d) 
is smaller than between (a) and (c) is the weaker phase 
transition in the case (d). The strength of the PhT is 
controlled by sq/sh, the ratio of entropy densities of the 
QGP and hadron phase at T c which is ~ 15 for the case 
(c) and ~ 3 for (d). If the ratio in the case (c) is artifi- 
cially forced to that in (d) by decreasing the number of 
degrees of freedom in the QGP, the changes in the pion 
spectra from (a) to (c) and (b) to (d) become approxi- 
mately the same. We conclude, that for the realistic EoS 
with HRG and QGP, the change in the slope of thermal 
pions, is caused mainly by heavy hadrons and resonances 
in the EoS and to much lesser degree by the phase tran- 
sition. 

The values of the average transverse momentum (pr) 
and multiplicities of positive thermal pions are collected 
to the Table |. In addition, the column (e) shows the 
same quantities for HRG+QGP when also the pions from 
the resonance decays (RD) are included. The results 
show that (pt)'s of thermal pions increase ~ 10 % when 
Tdcc drops from 160 to 100 MeV in all the cases (a)-(d). 
However, the change is ~ 21 % in (e), where 66 % of 
the pions come from the decays at Td cc = 160 MeV (cf. 
Table Q). In this case the low-pr part of the pion spec- 
trum, that gives the dominant contribution to multiplic- 
ity and (pt), is filled by decay pions from slowly moving 
resonances. As Td ec decreases the fraction of resonances 
drops increasing the growth of the average transverse mo- 
mentum. The same behaviour is seen in the case (b), if 



decay pions are included. The last three columns of Ta- 
ble | show the total multiplicity dN/dy of all particles 
and their total transverse energies dEx/dy and dEr /drj, 
where E T = E sin 9 || . 




p T [GeV] p T [GeV] 



FIG. 2: Transverse spectra of charged hadrons in Au+Au 
collisions at ^J~s — 130 AGeV with a 5 % centrality cut. The 
data is from PHENIX [jlij and the curves are our results for 
Tdcc = 120, 150 and 160 MeV. 

Comparison with data. Next we demonstrate that (i) 
using the initial conditions from the pQCD + saturation 
model and (ii) the EoS HRG+QGP with large number 
of hadron states in the low temperature phase, it is pos- 
sible to find a single decoupling temperature which re- 
produces both the measured abundances of pions, kaons 
and (anti)protons and the shapes of their spectra. 

Transverse spectra of pions, kaons and (anti)protons, 
calculated at T^ cc = 120, 150 and 160 MeV, are shown 
in Fig. 2 with the data from PHENIX Spectra of 

negative (positive) particles are in the left (right) panel. 
We find the overall agreement between the data and 
the calculations at Td cc = 150 . . . 160 MeV very satis- 
factory, considering that we have not adjusted parame- 
ters other than the decoupling temperature. A ~ 30 % 
feed-down from hyperons is included in the pip) data 
but not in the calculation. Especially, we would like to 
draw attention to the Td cc dependence of the slopes: at 
Tdcc ~ 150 . . . 160 MeV the slopes agree with the slopes 
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TABLE II: The multiplicity dN/dy, the fraction F of ther- 
mal particles in the multiplicity, and the average transverse 
momentum (px) in GeV for positive pions, kaons and protons 
calculated at three decoupling temperatures. The last column 
is the antiproton-to-proton ratio p/p = (dN p / dy) / (dN p / dy) . 



of the data but at 120 MeV all spectra are too shallow. 

In Table [n] we list the multiplicities and the thermal 
fractions F for positive pions, kaons and protons, also 
at T doc = 160, 150 and 120 MeV. Our results on the 
multiplicities agree within the experimental errors with 
the measured values |fL5f . For pions the thermal frac- 
tion F increases strongly because the number of heavier 
resonances drops fast. This stronger suppression of heav- 
ier states at lower temperature is the sole source for the 
relative change between the multiplicities of pions and 
heavier particles; it does not depend on the flow at all. 

The widening of spectra with decreasing Td oc is also 
seen m Table || as an increase in the average trans- 
verse momenta. For pions and kaons our results at 
T dec = 150 .. . 160 MeV are ~ 10 % larger than the 
measured values and slightly outside the experimental 
error bars [jl5|. The dEx/dr] which we obtain with 
Td cc = 150 . . . 160 MeV agrees within 8 % with the latest 
PHENIX result (jl6|, when the difference in the defini- 
tions of Et is taken into account. 

PHENIX has also reported p (p) yields corrected for 
hyperon feed-down, dN/dy = 19.3 ± 0.6 (13.7 ± 0.7) 
@. Our results in Table H for T dcc = 150 MeV agree 
with the measurement. With Td cc = 150 (160) MeV we 
also obtain A/A = 0.70 (0.74), again consistent with the 
PHENIX result 0.75 ± 0.09 0. Note that the devia- 
tion of p/p and A/ A from unity is due to the net baryon 
number content of the initial matter at saturation. 

We have not performed a x 2 fit for finding the best 
value for Td ec since our aim is to see how well we can 
describe the main features of the data while keeping our 
framework as simple as possible: pQCD + saturation to 
calculate the initial state, EoS to govern the expansion, 
and Tdec to specify the freeze-out. Also all the input 
parameters with the exception of Td GC are based on our 
previous studies j§ §, |l|l and we have not tried to tune 
them here. 

Summary. Using a hydrodynamic framework, we have 
studied the effect of the EoS and the decoupling temper- 
ature on the particle spectra and multiplicities. We have 
demonstrated an interesting and significant interdepen- 



dence between the hadron content of the EoS, the decou- 
pling temperature and the strength of the flow. Taking 
the initial conditions from the pQCD + saturation model 
|g, ||, we have calculated the transverse spectra and the 
multiplicities of pions, kaons and (anti)protons at the 
RHIC energy ^/s — 130 AGeV with a 5 % centrality cut 
in gold-gold collisions. The agreement between our re- 
sults with Tdcc = 150 . . . 160 MeV and the data measured 
by the PHENIX Collaboration Jl|, is very good both 
for the shape and the normalization of the spectra. 
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